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Abstract 

The evaluation of the absorption coefficients are important for particle emission caused 
by Hawking radiation. In the case of cosmological particle emission from the event 
horizon in De Sitter space, it is known that the scalar wave functions are solved in terms 
of Legendre functions. For fields with higher spin, the solution has been examined with 
low frequency approximation. We shows that the radial equations of the fields with 
spin 0, 1/2, 1 and 2 can be solved analytically in terms of the hypergeometric functions. 
We calculate the absorption probability using asymptotic expansion for high frequency 
limit. It turns out that the absorption coefficients are universal to all bosonic fields; They 
depend only on the angular momentum and not spin. In the case of spin 1/2 fermions, we 
can also find non-vanishing absorption probability in contrast to the previously known 
result. 
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Since the celebrated discovery of the evaporation of the black holes |1| much effort 
has been given to get analytic expression of the evaporation ratepj. The fact that the 
the corresponding field equations around black holes can not be expressed by the special 
functions prevent us from getting analytic expression for high-frequency mode. On the 
other hand, the Hawking radiation caused by the cosmological event horizon fU seems 
to be much easier to handle because it is known that the scalar field equation in static 
coordinates of De Sitter space can be solved by Legendre functions Q]. We can therefore 
convert the solution at the event horizon to our observing world by analytic continuation 
for all the frequency. Also the Bogoliubov coefficients of the mode functions from the 
global coordinates has also been obtained 0. As for the fields with spin, the solution was 
made only by low frequency approximation where they did not found any radiation for 
spinor fields. Our aim of this short letter is to show that the the radial equation for spin 
fields can also be solved by special functions. We obtain the absorption probability using 
high frequency approximation. It turns out that the absorption probability for bosonic 
fields are universal and do not depend on the spin. For spinor, we find non-vanishing 
result for the absorption probability in contrast to the result given in ref. M. 

De Sitter metric in static coordinates takes the form; 

2 2 

ds 2 = [1 - {%)}dt 2 - [1 - (^)rMr 2 - r 2 (dd 2 + sin# 2 V), (1) 
a 2 a 2 

where a is related to the cosmological constant by A = 3/a 2 . In the case of the type D 
space-time with vanishing electromagnetic fields, the equation of the first order variation 
of the scalar, spin 1/2 fermions, electromagnetic fields and gravity is universal and written 
by Teukolsky equations]!]] and separable. In our case, the angular component can be 
solved in terms of the spin-weighted spherical harmonics s Y{ n {9) given in Ref. |§. Then 
the radial equation is given in the variable z = r/a as|| 

z 2 {1 _ z 2 f{^) 2 R v - [2(p + l)z 3 - 2(s + 1)(1 - z 2 )z](l - z 2 )j-R p 
-{(1 - z 2 )[(l - s)(l + s + 1) + (s + p + l)(s + p + 2)z 2 } - (auozf + 2iauzp}R p 
= 0, (2) 
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where p takes the value ±s. In the case of the scalar fields s = p = 0, by choosing 
variable as y = z 2 , the singular point of the equation reduces to 0, 1 and oo with definite 
singularity. The equation can then be solved in hypergeometric functions^]. This choice 
of the variable is a natural choice because the metric depends only on z 2 . For s / 0, 
however, this choice does not work. There is another natural choice of variable. Near 
horizon, a useful coordinate is given by dz* = dz/(l — z 2 ) and therefore we choose a 
variable as 

» = = IT? (3) 

With this choice of the variable, Eq.(2) can be written as 

+ [<p + l)i - 2(s + 1 )T ±- - 2(. + P + D^l^ 

- w-^ + s + i)^-(- s+P+ i)(- s+P+ 2)^ + l -^i! + ^] 

+ [-(/ - s)(l + s + 1) - iacup] +(s+p+ l)(s+p + 2) 1 }R P = 0.(4) 

y{y -i) +i) 

Now the equation have four defenite singularity at points 0, —1, +1 and oo. It turns out 
that the singularity at y — — 1 can be factorized as 

Rp = y- p -^(i - y) l - s (i + y) s+p+l U (5) 

and f p satisfies the hypergeometric equation 

y{l - y)f; + [l-p-iaw-{2l + 3-p- iau)y}f p - (/ + 1 - p){l + 1 - iaw)f p = 0, (6) 

which has two independent solutions. Out of two solutions, we take the outgoing solution 
at the future horizon y ~ 0(z ~ 1) which behaves as exp(2p + iau)z* where z* = 
— (logy)/2. The solution is 

Rp = y-^~ p (i - y) l -\y + 1) s+p+1 F(a, b- c- y), (7) 
where F(a, b; c; y) is the hypergeometric function with 

a — I + 1 — p, b — I + 1 — iau, c = 1 — p — iau. (8) 



By using the recursion relation of the hypergeometric functions, we can derive the 
recursion relations which relate the solutions with different spins; 

^b--(l-y) 2s+1 (l + y)- 2s -^] 

= _( a + ia co)y-^(l - y) 2s+ \l + yr 2s - 3 R s+ i, (9) 

^[y^ +s+1 (i-y)(i + y)- 2s - 3 R s+ i} 

= - (/ '^ + g + 1) ^d - y)-\i + vY^Rs. do) 

s + mtu 

We can also obtain the operations which change the angular momentum. 

The behavior of this function near us (y ~ 1, z ~ 0) can be obtained by the analytic 
continuation. Near z — 0, the variable y can be written as y ~ e _2z , 1 — y ~ 2z. Since 
the parameters in the hypergeometric function satisfies a relation c—a — b = — (2/ + 1), the 
analytic continuation of the function produces logarithmic terms and the analysis becomes 
very complicated. Therefore, we will analytically continue the angular momentum to 
avoid such complication. (Strictly speaking, we cannot justify this procedure, although 
this technique is also used in the case of Black hole@.) Around z — 0, R p behaves as 



R P 



T(l-p-iauj)T(-(2l + l)) , , 2z 

JJ -F{l + l-p,l + l-iauj;2l + 2; 



T(-l -iau)T(-l-p) v " ' 'l + z' 

r(l - p - ia^)V(2l + 1) 2z_ rP8+1)F( _ z _ pj _ z _ iaw . _ 2/ . Jf_ )](11) 



r(/-j9 + i)r(/ + i-mcu) v i + z / v ™ 'i + z' 

We are going to consider the high frequency approximation au » 1. By using a formula 

F(a,b;c;x)~F(a;c;bx), (12) 

when 6 — > oo, and also considering the region where |aw2;| >> 1 and | z \« 1. Then we 
can use an asymptotic expansion of the confluent hypergeometric functions; 

FCicjxJ-pI^^'+S^x" (13) 
T(c — a) v x 7 r(a) 

and find the following asymptotic expansion of i? p ; 



We observe that i? p is not singular in the limit 2/ =integer. In order to calculate the 
absorption coefficients, we use a trick used in Refs. || |2], |5| . When we write the asymptotic 
expansion of R s as 

Rs ~ YSfi-*"/^* + Y^e^/z, (15) 
the absorption probability T can be calculated as[f| 

i_r = L , m ™ |. (16) 



ys y- 

- 1 Out - 1 Out 



From Eq.(14), we find 



Y Q s ut = A s A p {iau) l+p [l-{-l) 2p {iau)~ {2l+1) 



r(Z + 1 — zau;) 
T(— / — iauj) 



1 (— / — sj 1 (— / — zaa;) 



where 



A. = (-1) ^'}"'"!°"' . (18) 
By using Eqs.(16) and (17), we obtain 

4A Z 



(1 + A ; ) 2 ' 

where 



/ n 2 



(19) 



A; :inte ger = + 7 \2> ( l > S for bosons), (20) 



n=l 
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1+1/2 (n - ^ 2 



A i:ha if_i n te g er = Jl ( X + ~i — \2~) ( l ^ s for fermions). (21) 

n=i y auj ) 

It is amazing to observe that the absorption probablility is independent of the values 
of spins and depends on spices (fermion or boson) and the total angular momentum. 
For spin case, our result in Eqs.(19) and (20) agrees with the one derived by Lohiya 
and PanchapakesanfU] and the absorption probability vanishes for 1 = 0. For s = 1,2 
cases, we have non-vanishig probabilities for all I and our formula are different from those 
by Lohiya and Panchapakesan[6]. In particular, we found non- vanishing probability for 
s = 1/2 in Eqs.(19) and (21) which is in contrast to the result in Ref. |]] who argued that 
the absorption probability is zero in their approximation. 

For s = 0, the functions R p is the same as R~ p so that the trick used to derive T 
is not applicable. Lohiya and Panchapakesan|Q] derived T and showed that the result 
agreed with ours in Eqs.(19) and (20). As for s = 1/2, i?_ p is expressed in terms of R p 
so that another method using currents may be useful. For this purpose, we can use the 
conserved current to obtain the absorption coefficients as follows. We define 

g p = y^/\l - yY s+1 \l + y)^ + ^R p (22) 

and then we find that 

W = 9 ~ p ^~Jy 9 ^ ~ 9 ^~Jy 9 ^ = constant i ( 23 ) 
which is the Wronskian. For the solution in Eq.(7) for spin 1/2 case, we have 

9i/2 = y-^-^(l-y) l+1 F(l + ^,l + l-iaw,-iau + ^,y), 
9-1/2 = y-^-i\l-y) l+1 F(l + p + l-iaiu,-iaw + ~-,y). (24) 

Since the functions gi/2 and (7-1/2 are not independent, we can rewrite the conservation 
of current in the following form; 

W = y~ 1/2 (l - y)-\9l/29il2 - ^ 2 \f i 9-i/29-i/2] = \ + iau, (25) 

(Jj UJ \~ A £ 
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where the value of W is fixed by evaluation near the horizon (y ~ 1). We evaluate W 
by using the outgoing flux at the region z « 1 and auz » 1 and can obtain the 
absorption coefficient. 

We have used the analytic continuation of the angular momentum. But we cannot 
justify the procedure of the analytic continuation whereas the justification can be achieved 
for the evaluation of the absorption coefficients of black holes [fTOfl . At present, we cannot 
obtain the asymptotic expansion without using the analytic continuation. We hope that 
more rigorous approach can be done. 
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